Nodal solitons and the nonlinear breaking of discrete symmetry 
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We present a new type of soliton solutions in nonlinear photonic systems with discrete point- 
symmetry. These solitons have their origin in a novel mechanism of breaking of discrete symmetry 
by the presence of nonlinearities. These so-called nodal solitons are characterized by nodal lines 
determined by the discrete symmetry of the system. Our physical realization of such a system is a 
2D nonlinear photonic crystal fiber owning Cq v symmetry. 



Symmetry is one of the most powerful and elegant con- 
cepts in physics. The use of group theory provides an 
extraordinary mathematical tool to classify solutions ac- 
cording to the symmetries of the physical system. In 
recent years, the increasing interest in physical systems 
owning 2D discrete symmetries, such as 2D nonlinear 
photonic crystals or Bose-Einstein condensates in 2D pe- 
riodic potentials, raises the question of utilizing group 
theory as an analysis tool. Certainly, this approach is a 
standard in solid state physics. Its use in the topic of 
photonic crystals is less extended and it has been tradi- 
tionally confined to the classification of linear modes pj. 
Its generalization to the nonlinear case has become of 
great interest after the recent experimental observation 
of fundamental and vortex solitons in optically-induced 
2D nonlinear photonic crystals 0- S Q| - In fact, attempts 
to apply group theory to the analysis of this type of solu- 
tions permitted the analytical prediction of the angular 
dependence of vortex solitons in 2D photonic crystals 0. 
Following this approach, in this paper we will use group 
theory as a general framework to analyze the role played 
by nonlinearities in the realization of discrete symme- 
try. We will see that discrete symmetry is realized dif- 
ferently when nonlinearities are present and that a new 
phenomenon of discrete-symmetry breaking occurs. The 
physical outcome of this general process is the generation 
of a new type of solitons with lesser symmetry than that 
of the original system. 

So that, we start by analyzing the general problem of 
finding stationary solutions -$(x, y, z) — <fi(x, y) exp i/3z— 
of a nonlinear operator of the form: 
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where Lq is a linear operator (depending on the trans- 
verse coordinates x t = (x,y) only) invariant under a 2D 
discrete point-symmetry group G and Lml is a nonlinear 
operator depending locally on the modulus of the 4> field. 
We are interested then in solving the following nonlinear 



eigenvalue problem: 

{LQ + L NL {\ ( j>\)) ( j> = 2 ( j>. 



(2) 



Since Lq is such that [Lq,G] = (i.e., it is invariant 
under the action of all the elements of the group G: 
g Log = Lq, \/g 6 G), all its eigenmodes have to lie on 
finite representations of the discrete group G [f| . What 
we try to determine next is the effect of the nonlinear 
term Lml in the symmetry properties of Eq.{2J solutions. 

A solution 4> s of Eq.(2) nas to satisfy the so-called self- 
consistency condition, namely, <p a has to appear as an 
eigenmode of the operator generated by itself, L ((/>„) = 
Lq + Lnl(\4> s \). From a symmetry point of view, the 
self-consistency condition implies that if (f> s belongs to 
some representation of a finite group G', then the entire 
operator L(<fi s ) has to be invariant under the same group, 
[L((j) s ), G'] = (otherwise, L(<fr s ) would not contain in its 
spectrum the representation where <^ s lies on). We call 
this property the group self-consistency condition. 

The simplest attempt to find solutions of Eq.(2J satis- 
fying the group self-consistency condition is trying func- 
tions that enjoy the full symmetry of the linear opera- 
tor; i.e., functions that are invariant under G. Functions 
belonging to the fundamental representation of G sat- 
isfy this property Q: <^>f und = g0f un d = ^ftmd, V.g € G. 
Group self-consistency is satisfied because q~ 1 Lnl9 — 

WKndl) = Wl&andl), ^5 € Gd.e., [L NL ,G] '= 

and, thus, [L((f>{ un< i),G] = 0. Solutions that satisfy this 
property are called fundamental solitons and they have 
been found in different systems of the type described by 
Eq.JSJl- A less obvious choice is the selection of func- 
tions belonging to higher-order representations of the 
same symmetry group G of the linear system. For 2D 
point-symmetry groups, these higher-order representa- 
tions can be either non-degenerated (one-dimensional) 
or doubly-degenerated (two-dimensional) 0. Vortex- 
antivortex solutions, appearing always as conjugated 
pairs (4> v , 4>%), belong to two-dimensional representations 
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Figure 1: Two nodal solitons for Z = 1 (A = 23 fim, a = 8 /im, 
7 = 0.006 and wavelength A = 1064 nm): (a)-(b) amplitude 
and phase, respectively, of the S nodal soliton; (c)-(d) ampli- 
tude and phase, respectively, of the A nodal soliton. Inset: 
schematic transverse representation of a PCF. 



of G [s|. We note that the modulus of a vortex solu- 
tion is a group invariant (this is a general property also 
fulfilled by the modulus of functions belonging to one- 
dimensional representations of G). Since \<f)„ \ = \4>v\, 
then g~ 1 L NL g = L NL (\(f> 9 v \) = L NL (\(j) v \), \fg G G and, 
consequently, [Lnl,G] = and [L((p v ),G] = 0. It is 
apparent that vortex solitons also fulfill the group self- 
consistency condition. 

Up to now, we have considered solutions belonging to 
representations of G, the symmetry group of the linear 
operator. For them, the self-induced nonlinear opera- 
tor enjoys the same symmetry as its linear counterpart: 
[Lnl(\<Ps\),G] = [Lo,G] = 0. However, we ask ourselves 
if solutions with different symmetry than that exhibited 
by G can also fulfill the group self-consistency condition. 
Let us assume a trial function belonging to a certain rep- 
resentation of a group G' such that G' ^ G. Since the 
modulus of the function is G'-invariant, \(f)f \ = \<f) s \, the 
nonlinear operator is also G'-invariant, [£jvz,(|</> s |), G'} = 
0. If G' > G, the total operator L(4> s ) cannot have the 
G' symmetry because the linear operator has less sym- 
metry. The linear part of L(<j) s ) breaks the G' symmetry 
of the nonlinear part and the group self-consistency con- 
dition cannot be satisfied: [L(0 s ),G'] ^ 0. Thus, we 
disregard functions with symmetry higher than G. This 
scenario changes if one considers functions with lesser 
symmetry than G; more specifically, functions belonging 
to representations of a subgroup G' of G (G' C G) . The 
difference now with respect to the previous case is the 
following: since the linear part is invariant under the G 
group, [Lq, G] = 0, it is also invariant under any of its 
subgroups. Thus, it is also true that [Lq,G'] = 0. Since, 



as before, \4>f \ = \4>s\, the nonlinear operator verifies 
[Lnl(\<Ps\), G'] = and, consequently, the total opera- 
tor is G'-invariant, [L(<j> s ),G'] = 0. The function then 
fulfills the group self-consistency condition for the sub- 
group G'. Therefore, this type of functions can also be 
solution of the eigenvalue Eq.J5J (note that solutions with 
no symmetry also verify the group self-consistency condi- 
tion since the identity transformation constitutes a sub- 
group of any group G) . We cannot guarantee that they 
are indeed solutions until we solve Eq.J3J with the con- 
straint 4> S D(G') (D(G') being a representation of G'), 
since the trivial solution (0 = 0) is always valid. In this 
sense, group self-consistency is a necessary but not suf- 
ficient condition. However, if these solutions are found, 
they would be characterized by the symmetry breaking 
pattern G — > G' C G. In other words, we would say 
that, for those solutions, nonlinearity breaks the discrete 
symmetry of the linear system. 

Our specific physical system is a triangular photonic 
crystal fiber (PCF), although similar results can be ex- 
pected in other 2D photonic crystals. We study the prop- 
agation of the electric component of a monochromatic 
electromagnetic field (at fixed polarization: E = 0u, 
|u| = 1). PCF's are thin silica fibers possessing a reg- 
ular array of holes extending the entire fiber length and 
characterized by the hole radius a and the spatial pe- 
riod A of the photonic crystal cladding (see inset in 
Fig. When silica nonlinearity is not neglected, a 

PCF is a particular case of a 2D nonlinear photonic 
crystal with a defect (where guidance occurs). In this 
case, Lq = Vj + kgno(x,y), where V* is the transverse 
gradient operator, is the vacuum wave number, and 
no is the refractive-index profile function (no = nr s iu ca \ 
in silica and no = 1 in air). The nonlinear term is 
Lnl — fcg7A(x, y)\4>\ 2 , where A is the distribution func- 
tion of nonlinear material (A = 1 in silica and A = in 
air) and 7 is a dimensionless nonlinear coupling constant, 

7 = 3 X( s ii ica )-P/ (2eocn(siiica)-4o) (P is the total power and 
Aq is an area parameter: A® = 7r(A/2) 2 ). 

The symmetry group of a triangular PCF is Ce v , i-e., 
[I/o,C6„] = 0. This group is constituted by discrete n/3- 
rotations plus specular reflections with respect to 

the x and y axes: 9 -4 —6 and 8 — + 7r — 0, in polar co- 
ordinates. Solutions with the PCF C6i,-symmetry have 
been previously found in the form of fundamental and 
vortex solitons 0, Q] . We will focus now on new solu- 
tions belonging to the subgroup C<2 V of Cev, formed by 
Rn,R x and R y . Thus, we study the particular sym- 
metry breaking pattern G = Cq v — ► G' = C% v , We 
can explicitly construct functions belonging to the four 
non-degenerated representations of Ci v out of functions 
in the two-dimensional representations of Cq v . The lat- 
ter functions come in conjugated pairs (<j)i, </>*) (I — 
1,2), whose angular dependence is fixed by symmetry: 
4>i = r e <j>f(r, 9) expfi^f (r, &)], where 4> s (r, 0) is a scalar 
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function, characterized by 4> s (r, 9 + 7r/3) = 4> s (r,9) and 
4> s {r,-9) = 4> s (t,tt - 9) = 4> s {r,9), and (j) p {r,9) is a 
pseudoescalar function characterized by (fP(r,6 + 7r/3) = 
</> p (r,<9) and <^(r, -0) = c^(r, tt - 6) = -<^(r,0). Let us 
consider the two following linear combinations (I = 1,2): 
l/-\/2[0z ±<fii]). By writing the angular dependence of 
4>i, the new functions adopt the form: 

4(r, 9) = V2rtf(r, 9) cos [19 + (r, 6) + 6] , I = 1,2. 

(3) 

where 8 = 0,7r/2. The four different type of solu- 
tions given by Eq.© belong to the four different one- 
dimensional representations of the Civ group; i.e., in 
Hamermesh's notation: 0q € e ^2,^0 

^!'^w/2 e According to our previous general ar- 

gument, the 4>\ functions can be solutions of Eq.Q for 
a PCF. Therefore, we solve Eq.l(2j with the ansatz given 
by Eq. (j^J by means of the Fourier iterative method pre- 
viously used in Refs.Q and [j| to find fundamental and 
vortex soliton solutions in PCF's. This method ensures 
that the group representation is preserved in the iterative 
process. Starting from a seed function of the form 10, 
the method either finds the trivial solution or converges 
to a solution belonging to one of the representations of 

Civ 

Solutions of Eq.(|2jl of the form given by Eq.jSJ are 
indeed found. They are characterized by nodal lines 
determined by symmetry through the implicit equation 
cos [19 + (r, 6) + 6] = (I = 1, 2), For this reason, we 
call them nodal solitons. The solution with S = cor- 
responds to the symmetric (S) Cg v combination <j) l s = 
l/V2((/)i + (f)*) and that with 5 — n/2 to the antisym- 
metric (4) one 4> l A = i/V2(<j)i -(f)*) (1 = 1,2). Note that 
the (j>i function is not a solution of the Eq.(|2l because 
the superposition principle does not hold. This function 
can only be approximated by a vortex solution in the lin- 
ear regime (7 « 0). Although the complete structure of 
nodal lines could be rather intrincate, S and A nodal soli- 
tons are characterized by principal nodal lines: a single 
principal line for I — 1 solitons and two orthogonal prin- 
cipal lines for 1 — 2. In our simulations, we have found 
this four different types of nodal solitons. Nevertheless, 
we will show here results corresponding to the S and A 
nodal soliton solutions with I = 1 only. In Fig. ^ we 
show the amplitude and phase of S and A nodal solitons, 
respectively. As predicted by the nodal line condition, 
the 1 = 1, S nodal soliton presents a single vertical nodal 
line, whereas for the A soliton this line is horizontal. 

The role played by symmetries can be clearly envisaged 
by analyzing the diagram n so ivs. 7, where n so i — f3/ko 
for different soliton solutions. In Fig. |2we represent the 
curves corresponding to the S and A nodal solitons as 
well as the curves for vortex solitons reported in Ref.[5]. 
In the linear case (7 = 0), the superposition principle 
holds and, therefore, the symmetric and antisymmetric 
combinations of the linear modes 4>\ and (the linear 
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Figure 2: Effective index of a soliton solution, n BO \ vs the 
nonlinear coupling 7 for symmetric (dotted line) and anti- 
symmetric (dashed line) solitons and vortex and antivortex 
solitons with I = 1 (solid line). 
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Figure 3: Evolution of a diagonal perturbation of a S nodal 
soliton (\<f)'(x, 0, z)| 2 ) showing asymptotic stability. 

modes of vortex- type with I = 1) are degenerated so- 
lutions. In the linear case, all of them (S, A, (pi and 
modes) belong to the same representation of the Cq v 
group (I = 1, or E2 in Hamermesh's notation) and, for 
this reason, they all have the same effective index. The 
presence of the nonlinearity changes this scenario. It 
provides different options for the realization of the dis- 
crete symmetry. Vortex solitons realize the discrete sym- 
metry of the linear system, [X(|</>i|),C6i;] = 0, whereas 
nodal solitons break this symmetry into its C2 V subgroup, 
[Xd^s^l), Civ] = 0. Consequently, their corresponding 
eigenvalues are different since they are no longer related 
by the original symmetry that they enjoyed in the lin- 
ear case (7 = 0). Moreover, group theory predicts that 
the vortex and antivortex solutions must have the same 
effective index as they provide the same total operator 
L(|(/>i|) and they belong to the same representation (E2) 
of it. This is not the case for nodal solitons. 4>s and 
4>a are not in the same representation of d v ■ This fact 
immediately implies that its corresponding eigenvalues 
must be different. Curves in Fig. explicitly manifest 
this feature. For small nonlinearities, all curves appear 
as nearly-degenerated, but as the value of 7 increases, a 
growing gap between nodal solitons and vortices occurs. 
Even larger values of 7 permit to manifest the existing 
gap between the S and A nodal solitons. 
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In order to check the stability of nodal solitons we need 
to solve the evolution Eq. after perturbing the solu- 
tion: <t>'s a — + The stability analysis follows 
that of vortex solitons in Ref. 0. In there we intro- 
duced the concepts of diagonal and non-diagonal pertur- 
bations. For one-dimensional representations, a diagonal 
perturbation is defined as that that preserves the repre- 
sentation in which the solution lies on. In the present 
case — in which 4>s and 4>a belongs to the B\ and Bi 
one-dimensional representations of Civ, respectively — , 
this definition implies that <f>' s € B\ and <p' A € Bi. Ex- 
plicit examples of such perturbations are scaled solutions: 
$S A = (l + e )0S,A, e ^ 0. Evolution yields numerical ev- 
idence that nodal solitons are stable under such pertur- 
bations, as shown in Fig in which a diagonal perturba- 
tion (scaled solution) is applied. However, non-diagonal 
perturbations, taking the perturbed solution out of its 
original representation, provide instabilities. These in- 
stabilities are of the oscillatory type, as shown in Fig. 
Q] and they can be understood as a simultaneous os- 
cillation among modes belonging to all the representa- 
tions of Ci v . This instability pattern, however, shows no 
trace of pseudo-soliton collapse nor of transverse ejection 
of pseudo-solitons typical of Kerr nonlinearities in ho- 
mogeneous media treated in the paraxial approximation. 
This particular behavior was first observed in vortices in 
PCF's 0. Since the self-focussing instability seems to 
be rooted in the paraxial approximation a plausible 
explanation of its absence is the non-paraxial nature of 
evolution in this case. Absence of ejection can be qualita- 
tive understood by the inhibition of transverse radiation 
induced by the photonic crystal cladding. 

An interesting interpretation of nodal solitons is as 
interacting pseudo-solitons. It can be proven that a S 
nodal soliton can be written as <ps = 4>o(x+xo, y)~4>a{x— 
xo,y), </>o being a localized function in the fundamental 
representation of Cq v . In the case that 4>q is a sufficiently 
localized function (large nonlinear coupling 7 or strong 
lattice index contrast) , 4>o can be approximated by a fun- 
damental soliton solution. Then, a nodal soliton can be 
envisaged as a pair of weakly interacting pseudo-solitons. 
Like in an homogeneous medium this interaction is re- 
pulsive Q. In our case, nor 7 nor the index contrast are 
necessarely large, consequently, the soliton-soliton inter- 
action cannot longer be considered weak. However, the 
group theory antisymmetric decomposition in terms of lo- 
calized solutions in the fundamental representation of C§ v 
remains valid. In this way, the concept of nodal soliton 
generalizes the idea of interacting pseudo-solitons into a 
regime of strong particle coupling (intermediate 7 and 
index contrast). In the weak soliton- interaction regime, 
new soliton solutions have been recently found that can 
also be interpreted in the context of group theory re- 
ported in this paper (lflj . 
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Figure 4: Four snap-shots of the evolution of a non-diagonal 
perturbation of a S nodal soliton. In this case, an oscillatory 
instability occurs. 
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